Polarons and Molecules in a Two-Dimensional Fermi Gas 
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We study an impurity atom in a two-dimensional Fermi gas using variational wave functions for (i) 
an impurity dressed by particle-hole excitations (polaron) and (ii) a dimer consisting of the impurity 
and a majority atom. In contrast to three dimensions, where similar calculations predict a sharp 
transition to a dimer state with increasing interspecies attraction, we show that the polaron ansatz 
always gives a lower energy. However, the exact solution for a heavy impurity reveals that both 
a two-body bound state and distortions of the Fermi sea are crucial. This reflects the importance 
of particle-hole pairs in lower dimensions and makes simple variational calculations unreliable. We 
show that the energy of an impurity gives important information about its dressing cloud, for which 
both ansatze give inaccurate results. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.-d 



Since interactions between atoms can be tuned to 
essentially any value, cold atomic gases provide a 
unique opportunity for studying experimentally many- 
body physics in regimes that cannot be realized in other 
systems. Recently, much attention has been given to the 
problem of a Fermi gas with a low concentration of a 
second species, a so-called highly imbalanced gas (see, 
e.g., [lH3|). One fundamental problem is the nature of 
the ground state of a single impurity atom in a Fermi 
gas. For weak interspecies attraction, the ground-state 
energy is well described in terms of a state with an impu- 
rity atom dressed by a single particle-hole excitation of 
the Fermi sea, often referred to as a "polaron" |4|,|5|], while 
for strong attraction, a state based on a molecular pic- 
ture gives a lower energy 0, 0] . The transition between 
the two states is predicted to be sharp d, Q ■ 

It is natural to ask whether this picture persists in 
lower dimensions. This is of theoretical interest, since on 
general grounds one would expect quantum fluctuations, 
in this case the creation of many particle-hole pairs, to 
play an important role. In addition, the problem is on 
the verge of being investigated experimentally with the 
use of optical lattices In one dimension, a polaronic 
description gives qualitative agreement with known exact 
results [loj ], In this paper, we consider the case of two 
dimensions. We perform simple variational calculations 
based on the polaron and molecule pictures, and these 
predict that the polaronic state has the lower energy for 
all interaction strengths, in marked contrast to what hap- 
pens in three dimensions. For an infinitely massive impu- 
rity, the problem may be solved exactly, and the results 
show that the actual ground state incorporates aspects 
of both pictures: A two-body bound state is present for 
all coupling strengths, in addition to distortions of the 
continuum states. We show that the energy of an impu- 
rity gives important information about correlations in its 
vicinity and about mutual interactions between impuri- 



ties at nonzero density. 

Model We consider a uniform two-dimensional (2D) 
Fermi gas of atoms of species a, to which is added a 
single impurity atom of species b. The two species may 
be either different hyperfine states of the same element or 
different atomic species, in which case 6 may be bosonic 
or fermionic. The 2D confinement may be realized by a 
very tight trapping potential in the transverse direction 
Q and weak longitudinal trapping. For densities low 
enough that only s-wave interactions are important, the 
Hamiltonian reads 
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since a-a interactions may be neglected because of 
the Pauli principle. The single-particle eigenstates are 
(x|k) = e lk x /-\/K, where we take the system to be en- 
closed in a 2D box of volume V = l? with periodic 
boundary conditions. The single-particle energies are 
e£ = k 2 /2m r7 , where m a is the mass of species a (we 
take h — 1 throughout); ak(&k) annihilates an a(b) atom 
in state |k). The interaction is modeled as v(q) = gi for 
particle momenta less than a cutoff value A and zero oth- 
erwise. The coupling gi and the cutoff may be eliminated 
in favor of the two-body binding energy > fll| , 
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with the reduced mass fi = (m" 1 + jti^ 1 ) -1 . For eb <C 
A 2 /2/x, none of the results depend on the cutoff, and we 
take A — > oo at the end of the calculation. 

Polaron We describe the state using the variational 
ansatz |j| 
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Figure 1: (color online) Energy E as a function of es for the 
polaron ansatz (for the q — approximation, see text) and 
the dimer ansatz. (For es < 2ef, the dimer has momentum 
PM 7^ 0; the thin straight line indicates the dimer energy at 
Pm = 0.) 



which describes an impurity atom b (here at zero momen- 
tum) dressed by a cloud of particle- hole excitations of the 
ground state \N) of TV a-atoms. This state is an expan- 
sion of the true state in terms of numbers of particle-hole 
excitations, truncated at first order. Minimizing the en- 
ergy functional for this ansatz leads to equations for the 
energy of the state relative to the ground-state energy of 
N a-atoms, E = {H)m P — Eq N \ and the amplitudes (j)^: 



E 



q 



n q T^E + eW), 



^kq = 1 T q (£ + 4 g) ) 
<f>o V £-Ae kq 



(3) 



Here the T- matrix T^(E 
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v 1 SikKAl 1 _ n k) i E - Ae kq ) plays the role of an 

effective interaction, with Aekq = ^ — eq + e qlk> an( ^ 
n q = ®{kp — q) is the Fermi function, kp being the Fermi 
momentum. In the thermodynamic limit (keeping the 
majority density n a = N/V = k F /4ir fixed), this be- 
comes an integral equation, whose root E may be found 
numerically. 

For now, we will focus on the equal-mass case, m a = 
m = 2fi. The energy E(€b) obtained numerically is 
shown in Fig. Q] In the weak- coupling limit e# <C 
€p = k F /2m, the energy E ~ — 2ep/ h\{2ep/es) de- 
creases very rapidly with increasing es- This result 
can be interpreted as a mean-field shift E ~ n a g, 
where the density-dependent effective interaction g = 
— Ait/m ln(n a 47rS 2 /mes) is the leading term in T q for 
E,£b 0. In the strong- coupling limit 6b 3> £f, 
we find E ~ — 6b — ?]£f, with 77 rj 0.26. Some an- 
alytical insight into these results may be obtained by 
approximating Aekq = (k 2 — k • q)/m by k 2 /m. as 
has proven qualitatively correct in the 3D case [12j . 
This "q = 0" approximation gives the equation E as 
n a To(E) = — cb exp(2eir/|-E'|) + 2eF, whose solution is 
plotted for comparison in Fig. [TJ It gives the correct 
weak-coupling limit; in the strong-coupling limit it yields 



E ~ — es — 2ep/eB, which misses the term 0(€f). We 
have also studied the case of arbitrary mi, and find that 
the mass dependence of the polaron energy is weak if 
energies are expressed in units of k F /2fi. 

At first sight, it is somewhat surprising that for strong 
coupling the polaronic ansatz leads to an energy close to 
that of a molecule in vacuum. In this limit, the overlap 
probability of >Fp with the non-interacting ground state 
is given by Z = \4>o\ 2 ~ 2eF/es, which tends to zero. 
Thus the state is comprised mainly of holes and an (ab) 
dimer, as may be seen from Eq. ([5]), and the coefficients 
</"kq [Eq. ©] reduce to the wave function for a molecule in 
vacuo plus a hole. The leading contribution to the energy 
is thus the energy of a molecule in free space, since the 
hole has an energy of at most ep. 

Molecule In three dimensions, it has been demon- 
strated that a simple variational wave function based on 
a molecular picture gives a lower energy than the po- 
laron ansatz ([2]) for sufficiently strong attraction 0|, and 
we investigate whether this happens in two dimensions. 
In its simplest form, such an ab dimer with zero total 
momentum may be modeled using the trial state Q 
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which corresponds to a correlated ab pair in states with 
momentum greater than kp and a Fermi sea with N — 1 
noninteracting atoms. Minimizing the energy functional 
leads to the following equation for the energy Em (again 
relative to that of the iV-particle Fermi sea), 
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whose solution is 
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and ipk oc 0(fc— kp)/{e^— Em)- This result is simple be- 
cause in 2D the density of states is independent of energy: 
The energy of the dimer is shifted by the kinetic energy 
of the lowest state not Pauli-blocked. In this approxi- 
mation, the bare zero-momentum dimer state is always 
energetically less favorable than the polaronic solution, 
and there is no sharp transition (Fig. [TJ , in contrast to 
the 3D case. 

The question arises whether dimers with nonzero mo- 
mentum have a lower energy. The extension of (U) 
to describe a dimer with momentum p is ^m ~ 
S|k|>fc F a kkp-kl-^' — 1)) which leads to an equation 
similar to Eq. The solution is 
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where es = h 2 /2pLa\ defines the 2D scattering length a^- 

k 2 

For €b > 2^7 i^S the dimer energy is minimal for p = 0. 

k 2 

For €b < ^ ■ though, a dimer at p = 7^ is favor- 
able, because the kinetic-energy increase is outweighed by 
the reduced Pauli-blocking shift j2Cj|. For concreteness, 
consider the equal-mass case: The dimer momentum 
Pm = 2yfcp02 — l/ a 2 changes smoothly from pm = (at 
6b = 2ep) to — > &f as — > ep/2, where the energy 
reaches the continuum threshold Em = (Fig. [1]). In 
that regime, Em = + 2 v / Fb(v / 2Tf — \Je~B~). Thus, the 
bare dimer is energetically unfavorable compared with 
the polaronic ansatz for any coupling strength. To ob- 
tain more insight, we now examine an exact solution. 

Infinite-mass limit We consider the case where the 
impurity atom is very massive, m& 3> m a . In addition to 
being exactly soluble by use of Fund's theorem this is 
also relevant to experiments with atomic mixtures (e.g., 
Li and Yb [13 )• The massive impurity may be treated 
as a static defect, which has two effects on the Fermi gas: 
It creates a two-body bound state with energy — e_e, and 
it also shifts the energy levels in the continuum. In any 
number of dimensions, the energy shift of an s-wave level 
in the continuum is given by —k8o{k) / m a R, where So(k) 
is the s-wave phase shift and R the radius of the sphere 
containing the Fermi gas. The density of s-wave levels 
per unit interval in k is n/R. The total energy change 
when an impurity is added to the iV-atom system is 

E = -e B — I F dkk5 (k). (7) 

m a Tr J 

In 2D, the scattering phase shift is given by cot 8o(k) = 
ln(/c 2 /2m a eB)/7r in the zero-range limit [111 ). 

For large e B , AE = E + e B ^ -e F + ln(e £ / /eF) . The 
leading term corresponds to the fact that the phase shift 
is close to 7r, and consequently each of the continuum 
s-levels has become close to the next lower one in the 
absence of the impurity, thereby lowering the energy by 
ep. The logarithmic term may be thought of as an ef- 
fective atom-dimer repulsion. This is in stark contrast 
to the result AEm = k 2 pj2m — > + (mj — > 00) ob- 
tained from the dimer ansatz. Moreover, for the po- 
laron in the heavy-impurity limit, we find numerically 
AE w — O.lAep + 0(n 2 ), which is much smaller in mag- 
nitude than the exact result, — ep. 

The exact result shows that the ground state has a 
bound two-body state and that levels in the continuum 
are modified. It therefore incorporates aspects of both 
the polaron and the dimer pictures. The missing ingre- 
dient in the trial states we have used is components with 
higher number of particle-hole pairs: This is treated only 
to first (zeroth) order in the polaron (dimer) ansatz. In 
the language of perturbation theory, the T-matrix em- 
ployed in the polaron ansatz is calculated in the ladder 
approximation, and thus contains only particle-particle 
and hole-hole scattering. However, it is known that 



the T-matrix for scattering of a fermion from a mas- 
sive impurity is essentially independent of the presence of 
the medium, because Pauli blocking of ladder diagrams 
is compensated by impurity-hole scattering (l5|. Our 
calculations show that the higher-order impurity-hole 
scattering processes change AE to leading order from 
— 0.14eF (for one particle- hole pair) to — ep in the exact 
result. 

However, the bare-dimer picture worked reasonably 
well in 3D, predicting a molecular transition in agreement 
with Monte-Carlo results. To understand this paradox, 
let us look at the role of dimensionality in the dimer prob- 
lem. Solving |5|) in D = 1, . . . , 3 (for equal masses) yields 
an energy of the form E M — —(b-^f+cd^f(^-) ( ~ D ^ 2 ^' 2 
as £b — > 00. The last term corresponds to an upshift 
of the dimer energy due to Pauli blocking of the states 
|q| < hp. In D = 3 the upshift vanishes in the limit 
e_B 3> 2tp: This is because the density of states g(e) oc \fk 
vanishes at low energies, so that the contribution from 
Pauli blocking for e < 2ep has a negligible weight for 
eb — > 00. The situation is dramatically different in lower 
dimensions. In D = 2, the density of states is constant 
and thus leads to an interaction-independent displace- 
ment of the vacuum energy by 2 x ep, recovering the 
total shift +ep. This also illuminates why bare dimers 
should be even less favorable in ID where the low- 
lying states have an even stronger weight g(e) oc \j\fi, 
leading to a diverging upshift for strong coupling. 

Dressing cloud of an impurity Important information 
about the structure of the dressing cloud of an impurity 
may be extracted from the results for the energy. As 
is done in the theory of dilute mixtures of helium iso- 
topes firl | . we define the quantity v — {dn a /dnb)^ a — 
~(d(ib/dn a )/(dfi a /dn a ), where n a is the chemical poten- 
tial of a er-atom. Physically, this is the number of a-atoms 
in the dressing cloud of an impurity. The requirement 
that fj, a be held fixed ensures that far from the impurity, 
the density of a-atoms is unchanged by addition of the 
impurity. For n <^ n a , this number can be deduced from 
the single-impurity energy, v = —dE/dep, which is plot- 
ted in Fig. O As expected, v tends to zero in the weak- 
coupling limit, v ~ 2/ ln(2eF/es), for m a = mj. For 
m& — > 00, we can infer that there is exactly one dressing 
atom as cb 00, v 1. This contrasts with the polaron 
ansatz, which for m — > 00 predicts vp — > i] w 0.14 (fol- 
lowing a peak near cb = 2ep), illustrating that the single- 
particle-hole picture highly underestimates the impurity 
dressing. For comparison, the bare-dimer ansatz predicts 
the unphysical result i>m — — 1 1 amounting to a deficit of 
atoms in the dressing cloud due to Pauli blocking. 

Nonzero impurity density An intriguing question con- 
cerns the behavior at nonzero impurity density: Do the 
"dressed" impurity atoms behave as fermions, bosons, or 
neither of them? For weak attraction, it is not implausi- 
ble that the dressed impurities have the same quantum 
statistics as bare ones. By contrast, for strong attrac- 
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Figure 2: (color online) Number of majority atoms v in the 
dressing cloud of an impurity as a function of ca- 



tion, one may expect the basic degrees of freedom to be 
best described in terms of ab dimers, which are bosons 
for fermionic impurities and vice versa. On the basis of 
simple arguments, we cannot arrive at a definite conclu- 
sion about the statistics obeyed by the elementary exci- 
tations; to do so. it would be necessary to investigate the 
importance of exchange processes in a system with two 
impurities (2l| . 

Let us consider the case when the quasiparticles are 
fermionic. This could apply for weakly interacting 
fermionic impurities, but also for bosonic ones if they 
form a tightly bound dimer with a majority atom. We 
now show that the single-impurity findings have implica- 
tions for the thermodynamic properties at nonzero con- 
centration rib/ria <C 1. The total energy density £ of 
such a Fermi liquid then reads £(rib) — £(0) + £o(nb) + 
E(eB)nb + ^fn-l, where £(0) is the majority energy, 
£o( n b) = 7Tn b/ rn t denotes the kinetic-energy density, 
with the effective mass modified by interactions, and 
the term n^i^es) gives the energy reduction due to bind- 
ing of independent quasiparticles. Even in the absence 
of direct interactions between 6-fermions, there is an in- 
duced interaction between them, mediated by the major- 
ity Fermi gas It turns out to be repulsive owing 
to the Pauli principle and is characterized by the Lan- 
dau parameter / = v 2 deF/dn a . Note that, since in 2D 
the density of states dn a /deF = "ia/27r is constant, / 
is nonzero for n a — > 0. We mention that for bosonic 
quasiparticles, the effective interaction follows in a simi- 
lar fashion [19j , the difference being that there is a direct 
s-wave interaction and that the induced interaction is 
attractive. However, how that influences the induced in- 
teractions depends nontrivially on the degeneracy of the 
bosons and is left for future studies. 

Finally, we mention that the properties of impurities 
can be probed using techniques similar to those in 3D. 
By exciting collective oscillations, the effective mass is 
accessible [3j. With increasing coupling, this tends to oo 
for the polaron and to M for the dimer ansatz. Another 
key tool is radio- frequency spectroscopy 0|, where a b 
atom is transferred from its initial hyperfine level to an 



empty one via a pulse of frequency u>. For a polaron, 
the transition rate T(uj) decomposes into a quasiparticle 
peak oc Z8{uj — \E\), indicating the polaron contribution, 
and an incoherent background ri nc (w), which increases as 
(w - \E\f/ 2 for < u - \E\ < e F and falls off as for 
uj — \E\ ^ tp , if final-state interactions are ignored. This 
contrasts with the dimer ansatz, which yields rjj(w) oc 
Q(lu — \E\)/u> 2 without any quasiparticle peak. 

In summary, using variational wave functions, we find 
no evidence for a sharp transition between the polaron 
and the molecular picture in two dimensions. Compari- 
son with the exact result for a heavy impurity shows that 
both ansatze lead to inaccurate results for the dressing 
cloud of the impurity in the strong-coupling limit. This 
reveals the key role of many particle-hole pairs, and it 
reflects the importance of quantum fluctuations in lower 
dimensions. We conclude that more work is needed to 
understand the nature of the ground state of an impu- 
rity with finite mass in a two-dimensional Fermi gas. 
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